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It  is  found  that  the  solid  solubility  of  Ce  in  Nd2_xCexCu04±(5  is  limited  up  to  x  =  0.2.  A  semicon¬ 
ductor  to  metallic  transition  is  observed  at  600  °C  in  d.c.  conductivity  data,  which  coincides  with 
a  transition  in  temperature-dependent  area-specific  resistance  (ASR).  Ndi.8Ceo.2Cu04±«5  is  thermody¬ 
namically  and  chemically  stable  against  gadolinia-doped  ceria  (GDC)  up  to  1200  °C.  On  the  other 
hand,  it  reacts  with  a  yttria-stabilized  zirconia  electrolyte  to  form  Nd2Zr207.  At  700  °C,  the  ASR  of  a 
Nd! .8 Ce0.2 Cu04±($ /GDC/Nd! .8 Ce0.2 Cu04±<5  cell  sintered  at  800  °C  is  0.13  ohm  cm2,  and  the  ASR  proportion¬ 
ally  improves  with  increase  in  the  sintering  temperature  of  the  electrochemical  cell.  The  improved  ASR 
and  electrochemical  performance  are  attributed  to  the  nanocrystalline  nature  of  the  cathode  material. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs),  which  convert  chemical  energy 
directly  into  electrical  energy,  have  been  viewed  as  promising  new 
power-generating  systems  and  true  multi-fuel  energy  devices  [1]. 
Conventional  SOFCs,  however,  exhibit  major  disadvantages,  such 
as  cell-cracking,  the  formation  of  non-conducting  compound(s) 
at  electrode/electrolyte  interfaces,  and  the  need  for  expensive 
non-metallic  interconnects.  The  drawbacks  are  mostly  due  to 
extremely  high  operating  temperatures.  Consequently,  interme¬ 
diate  temperature  solid  oxide  fuel  cells  (ITSOFCs)  have  gained 
increased  importance.  Nevertheless,  achieving  high  cell  perfor¬ 
mance  at  reduced  temperatures  has  posed  a  great  challenge  for 
engineers  and  scientists  engaged  in  SOFC  development  activities. 

In  the  past,  the  utility  of  the  cathode  material,  lanthanum  stron¬ 
tium  manganite  (LSM),  was  limited  due  to  its  chemical  reactivity 
with  yttria-stabilized-zirconia  (YSZ)  [2].  Furthermore,  it  has  been 
reported  that  the  performance  of  LSM  reduces  rapidly  on  lower¬ 
ing  of  the  operating  temperature  due  to  its  poor  ionic  conductivity 
at  such  temperatures.  Recently,  Lai_xSrxCoi_yFey03_(5  (LSCF)  com¬ 
pounds  have  been  proposed  as  alternative  cathode  materials  for 
ITSOFCs  [3,4].  The  utility  of  LSCFs  has,  however,  been  limited  due 
to  the  lack  of  chemical  stability  and  a  thermal  expansion  coefficient 
(TEC)  that  is  incompatible  with  YSZ. 
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A  potential  cathode  material  must  exhibit  electronic  as  well 
as  oxygen-ion  conductivity  to  increase  the  triple-phase  boundary 
(TPB)  to  enhance  the  efficiency  of  SOFCs.  In  addition,  the  material 
should  exhibit  good  catalytic  activity  and  compatible  mechanical 
properties.  Consequently,  mixed  ionic  and  electronic  conductors 
(MIECs)  have  attracted  a  great  deal  of  attention.  Many  mixed  oxides 
of  the  form  A2BO4  (A  =  rare  earth,  alkaline  earth;  B  =  transition 
metal)  that  crystallize  with  the  tetragonal  I<2NiF4-type  structure 
(space  group  I4/mmm )  have  received  renewed  interest.  Such  A2B04 
oxides  have  been  of  considerable  interest  due  to  adequate  oxy¬ 
gen  conduction,  strong  electrocatalytic  activity  towards  the  oxygen 
reduction  and  TECs  comparable  with  those  of  conventional  elec¬ 
trolytes  [5,6].  The  I<2NiF4-type  structure  has  an  outer  intergrowth 
of  alternating  perovskite  (AB03)  and  rock-salt  (AO)  layers  stacked 
along  the  tetragonal  c-axis.  Excess  oxygen  ions  can  be  incorporated 
in  the  form  of  interstitial  species  into  the  rock-salt  layer,  without 
the  compensation  by  cation  vacancies  observed  in  the  AB03  type 
materials  [7-10]. 

Significant  mixed  ionic  and  electronic  conductivities  with 
appreciable  electrocatalytic  activity  have  been  found  in  Ln2Ni04 
and  Ln2Cu04  systems  that  have  the  K2NiF4-type  structure  [11-16]. 
In  addition,  the  thermochemical  stability  of  K2NiF4-type  mangan- 
ites,  ferrites,  cuperates  and  nickelates  has  been  reported  to  be 
higher  than  those  of  the  corresponding  perovskite-type  oxides  [15]. 
The  Co-  and  Fe-containing  K2NiF4  compounds  have  recently  been 
studied  extensively  [17,1 8  ].  Although  the  total  conductivity  and  the 
TEC  ofLa2_xSrxCu04±(5  materials  are  fairly  satisfactory,  the  materi¬ 
als  exhibit  rather  high  area-specific  resistances  (ASRs)  [16,19].  On 
the  other  hand,  the  ASR  of  a  Ndi  7Sr0.3CuO4  cathode  on  a  samarium- 
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doped  ceria  (SDC)  electrolyte  has  been  found  to  be  satisfactory 
[11].  Preliminary  studies  of  Nd2_xCexCu04±(5  (0  <  x  <  0.2)  with  both 
Ceo.9Gdo.1O1.95  (GDC)  and  Lao.gSro.iGao.sMgo^O^  (LSGM)  elec¬ 
trolytes  have  reported  poor  performances  [10].  The  authors  have 
attributed  the  poor  electrochemical  performance  to  coarse  and 
uneven  particles. 

In  particular,  the  solid  solutions  Nd2-xCexCu04±($  (0<x<0.2) 
have  been  of  special  importance  due  to  extrinsically  created  elec¬ 
tron  charge  carriers,  after  doping  with  Ce.  Their  super-conducting 
transition  in  them  has  been  found  to  be  dependent  on  x.  The 
crystal  structure  of  Nd2_xCexCu04±5  is  P-type  tetragonal  and  has 
copper-oxygen  planes  of  Cu02.  In  these  planes,  the  copper  ions 
have  a  square  oxygen-ion  neighbourhood  (01  site),  and  a  sandwich 
is  formed  by  the  layers  of  oxygen  ions  (O2  sites)  interspersed  with 
neodymium-cerium  ions  [20].  Accordingly,  there  are  four  possible 
modifications  of  tetragonal  structures,  T,  T\  T*  and  O,  for  A2B04 
compounds  [21,22].  The  solid  solutions  of  Nd2-xCexCu04±5  have 
invariably  been  prepared  by  solid-state  reaction  at  very  high  tem¬ 
peratures  (>1200 °C)  for  long  periods  of  time  (>40h)  [11,15].  A 
close  scrutiny  of  the  literature  revealed  that  Nd2_xCexCu04±(5  has 
received  negligible  attention  for  cathode  applications  in  ITSOFCs. 

There  is  sufficient  scope  for  preparing  Nd2_xCexCu04±5 
(0<x<0.2)  oxides  through  chemical  routes  to  improve  the 
electrochemical  properties  for  ITSOFC  cathode  applications.  In 
this  work,  various  compositions  of  Nd2-xCexCu04±($  (0<x<0.25) 
are  prepared  by  means  of  the  acetate  combustion  technique. 
The  prepared  samples  are  characterized  using  X-ray  powder 
diffraction,  scanning  electron  microscopy  and  four-probe  d.c.  con¬ 
ductivity  measurements.  In  addition,  electrochemical  impedance 
spectroscopy  (EIS)  studies  are  carried  out  on  the  symmetric 
cells  Ndi.8Ce0.2Cu04±5/GDC/Ndi.8Ceo.2Cu04±(5  sintered  at  800  and 
1200°C. 


2.  Experimental 

The  initial  reagents,  namely  cerium  acetate,  neodymium  acetate 
and  copper  acetate,  were  procured  from  Aldrich  Chemicals,  USA 
and  have  a  purity  of  >99.9%.  All  the  reagents  were  dried  at  120°C 
for  6h  to  remove  any  traces  of  moisture  present.  The  reagents 
were  dissolved  in  double-distilled  deionised  water  separately  in  the 
required  stoichiometric  ratios.  All  the  solutions  were  then  mixed 
in  a  single  corning  flask.  The  homogeneous  aqueous  solution  was 
then  heated  at  1 00  °C  while  subjecting  it  to  a  constant  stirring  using 
a  rota-mantle  to  facilitate  evaporation  of  the  water.  The  result¬ 
ing  precipitate  was  then  charred  at  about  250  °C  by  introducing 
it  to  an  electric  muffle  furnace.  It  took  >12  h  to  complete  the  entire 
procedure. 

The  residue  obtained  was  then  pulverized  to  form  powders.  Pel¬ 
lets  of  diameter  and  thickness  9  and  1-2  mm,  respectively,  were 
obtained  by  uniaxially  compressing  the  ground  powder  at  a  pres¬ 
sure  of  3  tons  cm-2  with  a  Specac  (UK)  stainless-steel  die-punch 
and  hydraulic  press.  The  resulting  pellets  were  initially  calcined 
at  800  °C  for  4  h  in  an  electric  furnace.  Subsequently,  they  were 
crushed,  and  pellets  were  made  in  a  similar  manner  as  described 
above  before  final  sintering  at  1000°C  for  4h.  The  samples  were 
allowed  to  cool  to  room  temperature  in  the  furnace. 

All  the  prepared  samples  were  subjected  to  structural  charac¬ 
terization  by  X-ray  powder  diffraction  (XRD)  with  a  PANalytical 
X’pert  PRO  (Philips,  the  Netherlands)  instrument  that  employed  Cu 
Ka  radiation.  A  curved  graphite  crystal  was  used  as  a  monochro¬ 
mator.  The  X-ray  diffraction  measurements  were  carried  out  in  a  20 
range  from  1 0  to  80°  with  a  step  size  and  time  per  step  of 0.020°  and 
5  s,  respectively.  The  XRD  data  were  refined  with  X’pert  Highscore 
plus  software.  The  sintered  densities  of  all  the  samples  were  deter¬ 
mined  using  Archimedes’  principle.  The  microstructures  of  the 


sintered  samples  were  examined  with  a  JEOLJSM-6380A  scanning 
electron  microscope.  Microhardness  was  measured  by  the  Vickers 
indentation  technique  (HMV-2  microhardness  tester,  Shimadzu, 
Japan). 

A  thin  platinum  film  on  both  flat  surfaces  of  the  sintered  pel¬ 
let  was  obtained  by  d.c.  sputtering  and  resulted  in  good  ohmic 
contacts  for  d.c.  electrical  conductivity  measurements.  Prior  to  the 
conductivity  measurements,  the  sample  was  spring-loaded  in  a 
ceramic  cell  holder  (Amel,  Italy)  and  heated  to  700  °C  for  1  h  to 
homogenize  the  charge  carriers.  The  resistance  during  the  cooling 
cycle  was  measured  as  a  function  of  temperature  using  the  four- 
probe  method  with  a  computer-controlled  Keithley  6221  current 
source  and  a  2 182 A  nanovoltmeter  in  delta  mode.  The  tempera¬ 
ture  of  the  sample  during  the  measurement  was  controlled  with  an 
accuracy  of  ±1  °C  with  a  Eurotherm  221 6e  temperature  controller. 
The  tip  of  a  calibrated  thermocouple  was  kept  in  the  vicinity  of 
the  sample  to  measure  its  actual  temperature.  The  oxygen  partial 
pressure  around  the  cell  during  the  conductivity  measurements 
was  maintained  using  electronic  mass  flow  meters  and  con¬ 
trollers  (Teledyne-FIastings,  USA).  Nitrogen  was  used  as  the  carrier 
gas. 

For  electrochemical  investigations,  the  initial  slurry/ink  of  the 
Ndli8Ce0.2CuO4±(5  cathode  was  obtained  as  follows.  A  gram  of 
Ndi.8Ce0.2CuO4±(5  powder  was  mixed  with  3  wt.%  polyvinyl  buteral 
binder,  sodium  free  corn  oil,  and  ethyl  methyl  ketone.  The  mixture 
was  then  ball-milled  using  a  Pulverisette-6  (Fritsch  Germany)  for 

1  h  with  300  revolutions  per  minute.  Thirty-five  tungsten  carbide 
balls  of  diameter  1 0  mm  and  80-ml  bowl  of  the  same  material  were 
used.  The  GDC  (10  mol.  %  gadolinium-doped  ceria)  nano-powder 
(Aldrich,  USA)  was  pressed  in  a  similar  manner  to  obtain  pellets 
of  diameter  9  mm  and  thickness  1  -2  mm.  They  were  then  sintered 
at  1400  °C  for  6h.  A  sintered  density  of  about  96%  was  achieved. 
Both  flat  surfaces  of  the  GDC  pellet  were  then  roughened  with  #  60 
grid  paper  and  cleaned  with  acetone.  The  slurry/ink  of  each  cath¬ 
ode  material  was  then  spin-coated  at  3000 rpm  for  60s  on  both 
flat  surfaces  of  the  sintered  GDC  electrolyte  to  obtain  the  following 
symmetric  cell  configurations. 

Ndi.8Ceo.2Cu04±(5/GDC/Ndi.8Ceo.2Cu04±(5 

The  symmetric  cells  were  initially  baked  at  600  °C  for  2  h  to  burn 
out  the  organic  binders,  and  finally  sintered  at  800  and  1200  °C  for 

2  h.  The  cells  sintered  at  800  and  1 200  °C  are  abbreviated  as  cell-800 
and  cell-1200,  respectively. 

Electrochemical  impedance  measurements  on  the  sintered 
symmetric  cells  were  carried  out  using  a  computer-controlled 
Solartron  1255B  FRA  in  combination  with  a  Solartron  SI  1287 
electrochemical  interface.  Prior  to  the  impedance  measurement, 
the  spring-loaded  samples  were  heated  to  700  °C  (' T 1)  for  1  h  as 
shown  in  Fig.  1.  During  the  cooling  cycle,  a  dwell  time  of  30  min 
was  allowed  at  each  set  temperature  (I2,  T3,  T4,  ...)  to  ensure 
that  thermal  equilibrium  was  achieved  in  the  cells.  Ten  minutes 
prior  to  the  end  of  the  dwell  time  (ts  1  in  Fig.  1),  the  real  and 
the  imaginary  parts  of  the  impedance  were  measured  as  a  para¬ 
metric  function  of  frequency  and  temperature  in  the  range  of 
0.01-1  x  106  Hz  and  480-700  °C,  respectively,  with  a  signal  ampli¬ 
tude  of  50  mV.  The  concentration  of  02  gas  (Po2 )  was  monitored  as 
discussed  previously.  The  EIS  data  obtained  were  analyzed  using 
software  developed  by  Scribner  advanced  software  for  electro¬ 
chemical  research  and  development. 

To  evaluate  the  chemical  compatibility  of  the  Ndi>8Ce0.2CuO4±(5 
cathode  with  the  GDC  and  YSZ  electrolytes,  the  symmet¬ 
ric  cells  were  sintered  at  1200°C  for  48  h.  The  interface 
powder  was  then  subjected  to  XRD  to  identify  the  phases 
present. 
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Time  — ► 

Fig.  1.  Temperature  profile  during  EIS  measurements. 


3.  Results  and  discussion 

3.1.  X-ray  powder  diffraction 

Typical  X-ray  powder  diffraction  (XRD)  patterns  of 
Nd2_xCexCu04±($  with x  =  0, 0.20  and  0.25  are  shown  in  Fig.  2(a)-(c), 
respectively.  In  general,  all  the  diffracted  lines  were  broader  than 
the  usual  ones.  This  broadening  is  attributed  to  the  nanocrystalline 
nature  of  the  materials.  The  obtained  XRD  data  were  profile-fitted 
with  X’pert  Highscore  plus  software  before  indexing.  A  close 
look  at  Fig.  2(a)  reveals  that  all  the  characteristic  diffracted  lines 
match  closely  with  the  Joint  Committee  on  Powder  Diffraction 
Standard  (JCPDS)  data  (File  No.  01-084-2177)  that  correspond  to 
pure  Nd2Cu04.  Lines  corresponding  to  either  the  initial  reagents 
or  any  intermediate  compound(s)  are  not  present.  Similarly,  all 
diffracted  lines  appeared  appear  in  Fig.  2(b)  match  closely  the 
JCPDS  data  for  tetragonal  Nd2Cu04,  with  only  small  deviations. 
No  lines  due  to  initial  reagents  are  seen.  On  the  other  hand,  in 
addition  to  the  lines  associated  with  Nd2Cu04  (Fig.  2(c)),  a  few 
weak  characteristic  diffracted  lines  due  to  Ce02  are  present  in  the 
case  of  x  =  0.25.  Nevertheless,  the  absence  of  lines  corresponding 
to  either  pure  reagents  or  intermediate  phases  confirms  the  for¬ 
mation  of  single-phase  tetragonal  Nd2_xCexCu04±(5  solid  solutions 
up  to  x<0.2.  Furthermore,  the  cerium  precipitated  out  of  the 
Ndi.8Ce0.2CuO4±(5  matrix  for  x  =  0.25.  These  results  are  in  good 
agreement  with  the  results  of  Spinolo  et  al.  [23].  To  ascertain 
the  formation  of  a  solid  solution  and  its  consequence  on  the  host 
lattice  structure,  the  lattice  cell  constants  for  all  samples  under 
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Fig.  2.  X-ray  powder  diffraction  patterns  of  Nd2-xCexCu04±(5  with  (a)  x  =  0,  (b)  x  =  0.2 
and  (c)x  =  0.25. 


t 

>* 

CO 

c 
0) 
c 

1 

jQ 

< 

10  20  30  40  50  60  70  80 

28!  degree  - ► 


• 

(a) 

•   GDC 

°.  YSZ 

•  Jfi 

♦  -1 

l 

sldi  sCeo  2C11O4+5  (AC) 
Nd2Zr207 

**  *  * 

A  A.  _  jfc  *  A.  A 

(b) 

r&Lct-r-JgL 

<► 

r~* 

- 

■ 

Fig.  3.  X-ray  powder  diffraction  patterns  of  Ndi.8Ceo.2Cu04±($  with  (a)  GDC  and  (b) 
YSZ  sintered  at  1000  °C  for  45  h. 


study  were  determined  using  Unitcell  computer  software  [24]. 
The  obtained  lattice  cell  constants,  a  =  3.956  and  c=  12.208  A,  for 
pure  Nd2Cu04  are  in  good  agreement  with  the  reported  values 
of  a  =  3.944 A  and  c=  12.1 69 A  [25].  The  variation  of  lattice  cell 
constants  a  and  c  of  Nd2_xCexCu04±5  with  x  is  depicted  in  Fig.  3. 
It  can  be  seen  that  the  lattice  cell  constant  c  decreases  linearly 
up  to  x  =  0.2,  and  remains  unchanged  for  x>0.2.  These  results 
substantiate  the  solid  solubility  limit  up  to  x  =  0.2  discussed  above. 
On  the  other  hand,  there  is  a  negligible  change  in  the  lattice 
constant  a  with  increase  in  x.  The  Nd2Cu04  lattice  contraction 
along  the  c-axis  due  to  Ce  doping  has  been  observed  earlier  [20]. 
According  to  the  literature,  Ce  remains  tetravalent  and  occupies 
the  8-coordinated  Nd3+  site  in  a  T  type  Nd2Cu04  structure.  The 
ionic  radius  (rg  =  0.97  A)  of  8-coordinated  Ce4+  [25]  is  smaller  than 
that  of  the  8-coordinated  host  Nd+3  (rh  =  1.0lA)  [20].  Thus,  the 
partial  substitution  of  the  latter  for  the  former  leads  to  Nd2Cu04 
lattice  contraction.  Because  of  the  absence  of  apical  oxygen  and 
the  existence  of  a  single  Cu02  plane  in  the  unit  cell  of  the  V 
type  Nd2Cu04  crystal  structure,  the  Cu  ions  are  4-coordinated. 
Consequently,  the  coupling  between  the  Cu02  layers  along  the 
c-axis  is  expected  to  be  weaker  than  along  the  a-axis  and  thereby 
leads  to  an  anisotropic  lattice  contraction.  Voronin  et  al.  [20]  have 
given  a  detailed  justification  for  the  lattice  contraction  along  the 
c-axis  due  to  partial  replacement  of  Nd3+  for  Ce4+  in  Nd2Cu04. 

Table  1  compares  the  crystallite  size  of  all  compositions  under 
study,  determined  using  the  X’pert  Highscore  plus  software  based 
on  the  following  expression: 


0.9A 
P  cos  0B 


(1) 


where  t,  A,  and  are  the  thickness  of  the  crystallite,  X-ray  wave¬ 
length  and  Bragg’s  angle,  respectively.  Here,  ft  is  obtained  by: 

P2  =  pm2-Ps2  (2) 

where  and  /3S  are  the  measured  and  standard  full  width  at  half 
maxima  (FWHM),  respectively,  of  the  diffracted  line.  The  XRD  pat¬ 
tern,  obtained  by  running  the  experiment  on  a  standard  silicon 
sample  provided  by  PANalytical  (Netherlands),  is  used  to  estimate 


Table  1 

Comparison  of  crystallite  size,  lattice  strain,  sinter  density  and  microhardness  num¬ 
ber  of  Nd2-xCexCu04±5  with  x  =  0-0.25. 


Composition 

cs  (A) 

Density  (%) 

HV 

Nd2Cu04 

770 

88.56 

635 

Ndi  .95  Ceo.05  Cu04 

772 

84.00 

627 

Ndi.gCeo.i  Cu04 

1056 

80.57 

619 

Ndi.85Ceo.i5Cu04 

1842 

81.12 

622 

Ndi.sCeo.2Cu04 

1964 

80.92 

614 

Ndi.75Ceo.25Cu04 

1620 

89.33 

532 

Cs  -  crystallite  size;  HV  -  microhardness  number. 
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Fig.  4.  Variation  of  lattice  cell  constants,  a  and  c,  with  x  in  Nd2_xCexCu04±($. 

the  value  of  A  close  scrutiny  of  the  data  in  Table  1  suggests  that 
all  the  compositions  under  study  are  nanocrystalline  materials.  Fur¬ 
thermore,  the  crystallite  size  increases  with  an  increase  in  x  within 
the  solid  solubility  limit.  By  contrast,  the  grain  size  of  samples  with 
x  =  0.25  is  smaller  than  that  of  x  =  0.20.  The  precipitation  of  Ce02 
in  the  Ndi.8Ceo.2Cu04±(5  matrix  hinders  the  grain  growth,  thereby 
resulting  in  relatively  smaller  grains.  This  result  also  supports  the 
observed  solid  solubility  limit  up  to  x  =  0.2.  From  Table  1,  it  can  be 
seen  that  the  sintered  density  of  Nd2_xCexCu04±(5  decreases  withx. 
The  reduction  in  sinter  density  is  attributed  to  the  increased  crys¬ 
tallite  size  (Table  1).  A  close  inspection  of  Table  1  reveals  that  the 
microhardness  number  increases  commensurately  with  sintered 
density  of  Nd2_xCexCu04±(5.  By  contrast,  the  low  microhardness 
number  of  Ndli75Ceo.25Cu04±5,  in  spite  of  the  high  sinter  density, 
is  attributed  to  inferior  grain  glue/adhesion  resulting  from  the  bi¬ 
phase  mixture.  The  homo-junctions,  i.e.,  the  interfaces  between 
same  grains  (grain  boundary),  are  chemically  strong  and  thereby 
exhibit  better  adhesion  (glue)  between  the  grains.  On  the  other 
hand,  hetero-junctions  are  chemically  weakly  bound.  These  factors 
are  responsible  for  low  microhardness  of  the  bi-phase  composition, 
i.e.,  x  =  0.25. 

The  XRD  patterns  of  Ndi.8Ce0.2oCu04±(5 +  GDC  and 
Ndi.8Ce0.2oCu04±(5 +YSZ  sintered  at  1200°C  for  48  h  are  shown  in 
Fig.  4(a)  and  (b).  All  the  diffraction  lines  in  Fig.  4(a)  closely  match 
with  those  of  either  Ndi.8Ce0.2CuO4±5  or  GDC.  On  the  other  hand, 
in  addition  to  the  diffracted  lines  due  to  Ndi.8Ce0.2CuO4±($  and 
YSZ,  a  few  lines  corresponding  to  the  intermediate  compound 
Nd2Zr207  are  clearly  seen  in  Fig.  4(b).  These  results  suggest 
that  Ndi.8Ceo.2oCu04±5  exhibits  excellent  chemical  stability  with 
GDC  as  compared  with  YSZ.  Thus,  GDC  is  selected  for  further 
electrochemical  studies  of  the  Ndi  8Ce0.2CuO4±5  cathode. 

3.2.  Scanning  electron  microscopy 

Scanning  electron  microphotographs  (SEM)  of  Nd2_xCexCu04±(5 
with  x  =  0.05-0.2  and  0.25  are  displayed  in  Fig.  5(a)-(e),  respec¬ 
tively.  For  all  compositions  in  general,  90%  of  the  cleaved  surface 
consists  of  trans-granular  cleavages,  and  the  other  10%  consists  of 
voids  and  grain-boundaries.  A  close  look  at  the  figures  also  reveals 
that  the  grain  size  increases  with  increase  in  x  of  Nd2_xCexCu04±($. 
Furthermore,  the  grains  are  elongated  with  increased  cerium  con¬ 
tent. 


The  SEM  photograph  of  a  fractured  surface  across  the  elec¬ 
trode/electrolyte  interface  of  the  symmetric  cell-800  after  four 
thermal  cycles  is  depicted  in  Fig.  6.  As  seen  from  the  figure,  the 
electrode  and  electrolyte  form  intimate  contact  across  the  interface. 
Moreover,  electrode/electrolyte  interface  is  homogeneous.  Despite 
four  thermal  cycles  at  800-20  °C,  no  cracks  or  separations  between 
the  electrode  and  electrolyte  is  found  at  the  interface.  The  electrode 
layer  thickness  is  approximately  62.7  p,m.  A  close  look  at  Fig.  6(b) 
reveals  that  the  nano-crystallites  of  Ndi>8Ce0.2CuO4±($  have  agglom¬ 
erated  to  form  a  porous  electrode.  Furthermore,  such  agglomerates 
are  uniformly  distributed  in  such  a  way  that  the  electrode  structure 
becomes  highly  porous.  Sintering  of  a  symmetric  cell  at  1200°C 
facilitates  the  fusion  of  nanocrystalline  grains  of  Ndi.8Ce0.2CuO4±5, 
leading  to  grain  growth  (Fig.  6(c)).  By  contrast,  the  electrolyte  is 
nonporous,  i.e.,  highly  dense  (Fig.  6(a)). 

3.3.  d.c.  conductivity 

The  variation  of  d.c.  conductivity  with  temperature  for  all  the 
samples  under  study  is  shown  in  Fig.  7.  Nd2_xCexCu04±($  with 
x  =  0.1,  0.15  and  0.2  exhibits  a  positive  temperature  coefficient 
(PTC)  to  negative  temperature  coefficient  (NTC)  transition  at  about 
625  °C.  The  absence  of  the  transition  for  pure  Nd2Cu04  and  x  =  0.05 
is  in  line  with  other  reported  results  [21].  Below  this  transition 
temperature,  all  the  samples  obey  the  Arrhenius  law,  i.e., 

a  =  <70  exp  (2^j  »  (3) 

where  oq,  k,  T  and  Ea  are  the  pre-exponential  factor,  Boltzmann 
constant,  absolute  temperature  and  activation  energy,  respectively. 
The  variation  of  d.c.  conductivity  at  600  °C  and  activation  enthalpy 
as  a  parametric  function  ofx  in  Nd2_xCexCu04±($  is  shown  in  Fig.  8  . 
The  d.c.  conductivity  initially  increases  with  an  increase  in  Ce  con¬ 
tent,  and  attains  a  maximum  for  x  =  0.1 5  and  decreases  thereafter. 
The  activation  energy  follows  an  opposite  behavior  to  that  of  the 
conductivity  with  x.  The  minimum  activation  energy  composition 
coincides  with  the  composition  that  has  the  maximum  conductiv¬ 
ity. 

Intrinsically,  electrons  and  holes  can  be  produced  in  Nd2Cu04±(5 
lattice  by  the  following  two  point  defect  equilibria  [21  ]: 

00  V0**  +  2e'  +  lo2,  (4) 

lo2  O'l  +  2h*,  (5) 

coupled  with  intrinsic  ionization: 

nil  e'  +  h  (6) 

At  intermediate  Po2  and  temperature,  Nd2Cu04±5  is  almost  sto¬ 
ichiometric  in  oxygen,  and  shows  semiconducting  behavior  where 
the  carriers  are  electrons  and  holes  created  by  thermal  promo¬ 
tion  of  electrons  from  the  valence  band  to  the  conduction  band 
according  to  Eqs.  (5)  and  (6).  At  higher  temperatures,  however,  the 
majority  charge  carriers  are  electrons  created  by  the  equilibrium 
Eq.  (4)  in  the  whole  range  of  P02.  At  the  same  time,  there  may 
be  an  overlapping  of  valence  and  conduction  bands  leading  to  a 
metallic  behavior,  i.e.,  the  conductivity  decreases  with  an  increase 
in  temperature  (Fig.  7). 

The  partial  substitution  of  Ce4+  for  Nd3+  introduces  extrinsic 
defects  according  to  the  following  quasichemical  reaction: 

CuO  +  2Ce02  -*  Cue/  +  2CeNd#  +  400*  +  2e'  +  ^02  (7) 

The  electrons  thus  created  are  associated  with  the  Cu02  planes. 
The  concentration  of  such  extrinsically  created  electrons  in  Cu02 
planes  increases  with  an  increase  in  Ce  concentration.  The  increase 
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Fig.  5.  Scanning  electron  microphotographs  of  Nd2-xCexCu04±(5  with  (a)  x  =  0.05,  (b)  x  =  0.1 ,  (c)  x  =  0.1 5,  (d)  x  =  0.20,  and  (e)  x  =  0.25. 


in  electron  concentration  has  been  supported  by  Hall  coefficient 
measurements  by  Takagi  et  al.  [26].  Thus,  the  d.c.  conductivity  of 
Nd2Cu04±(5  in  the  present  study  increases  with  Ce  content.  Singh 
et  al.  [27]  have  observed  a  relatively  rapid  increase  in  the  con¬ 
tent  of  the  lattice  oxygen  of  Nd2_xCexCu04±(5  for  0.1 55  <  x  <  0.1 75. 
According  to  their  research,  there  exists  a  critical  range  of  unit  cell 
volume  and  oxygen  content  for  a  smooth  super-conducting  transi¬ 
tion.  Furthermore,  there  is  a  critical  unit  cell  volume  corresponding 
to  x  =  0.15.  Below  this  concentration  (x<0.15),  the  electrostatic 
attraction  amongst  the  lattice  is  strong  enough  to  limit  the  removal 
of  oxygen.  These  factors  may  be  responsible  for  the  maximum  con¬ 
ductivity  for  x  =  0.1 5  (Fig.  8).  The  variation  of  d.c.  conductivity  (at 
650  °C)  with  log  Po2  depicted  in  Fig.  9  suggests  that  the  electronic 
conductivity  decreases  with  an  increase  in  oxygen  partial  pres¬ 
sure.  The  very  small  d.c.  conductivity  that  is  dependent  on  oxygen 
partial  pressure  is  in  good  agreement  with  previous  works  [20]. 
Although  the  d.c.  conductivity  is  highest  for  Nd2_xCexCu04±(5  with 
x  =  0.1 5,  the  cathode  composition  with  x  =  0.20  is  selected  for  elec¬ 
trochemical  studies  with  GDC  as  the  electrolyte  because  the  solid 
solubility  of  the  GDC  electrolyte  is  up  to  x  =  0.2.  The  solid  solubil¬ 


ity  of  the  electrolyte  in  the  cathode  may  affect  the  electrochemical 
performance. 

3.4.  Electrochemical  impedance  spectroscopy 

Typical  complex  impedance  spectra  of  cell-800  at  660  and  480  °C 
are  depicted  in  Fig.  10(a)  and  (b),  respectively.  In  general,  con¬ 
voluted  impedance  spectra  are  observed.  The  impedance  spectra 
shown  in  Fig.  10(a),  analyzed  using  a  computer  program,  show 
a  clear  semicircular  arc  in  the  low-frequency  region.  The  incom¬ 
plete  semicircular  arc  in  the  high-frequency  region  is  attributed 
to  the  limitation  of  the  Solartron  1255B  FRA  for  high  frequencies 
(<1  MHz).  A  close  inspection  of  Fig.  10(b)  reveals  the  presence  of 
three  semicircular  arcs.  The  first  low-frequency  semicircular  arc 
is  due  to  the  passage  of  current  through  the  electrode  and  across 
the  electrode/electrolyte  interface  and  also  complex  electrode  pro¬ 
cesses.  The  second  semicircular  arc  in  the  mid-frequency  domain 
is  attributed  to  the  grain-boundary  contribution  of  the  GDC  elec¬ 
trolyte.  The  incomplete  high-frequency  semicircular  arc  is  due  to 
the  bulk  grain  of  GDC.  The  semicircular  arc  representing  elec- 
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Fig.  6.  Scanning  electron  microphotographs  of  (a)  electrode/electrolyte  interface  of  cell-800  and  cathode  surface  of  (b)  cell-800  and  (c)  cell-1200. 


trode  polarization  is  characteristic  of  the  reversible  reduction  of 
O2  gas  at  the  cathode,  and  subsequent  migration  of  O2-  across  the 
interface. 

Complex  impedance  spectra  of  cell-1 200  at  660, 680  and  700  °C 
are  given  in  Fig.  1 1 .  In  this  case,  two  discernible  semicircular  arcs  are 
present.  In  general,  the  real  axis  intercepts  of  all  semicircular  arcs 
increase  with  a  decrease  in  temperature  as  expected,  and  indicate 
increased  ASR  as  well  as  bulk  resistance.  The  electrical  equivalent 
circuit  shown  in  Fig.  12,  represents  the  impedance  data  for  all  sym¬ 
metric  cells  under  study.  The  constant  phase  element  (CPE)  in  each 
circuit  accommodates  the  depression  in  the  semicircular  arc  due 


to  the  distributed  relaxation  time.  Riead,  Select  and  Rcath  represent 
the  lead  resistance,  the  resistance  of  the  grain  and  grain  boundary 
of  electrolyte  and  the  resistance  of  the  cathode,  respectively.  Celect.  * 
and  Cdl  are  the  capacitances  due  to  the  grain,  grain  boundary  and 
double-layer  interface  capacitor,  respectively.  The  W  in  the  parallel 
circuit  representing  the  low-frequency  arc  is  due  to  the  Warburg 
impedance. 

The  electrolytic  bulk  conductivity,  crGDC,  at  atmospheric  pres¬ 
sure  is  obtained  using  Eq.  (8), 

aGDC  =  RS'  W 

where  t  and  S  are  the  thickness  and  the  surface  area  of  the  elec¬ 
trolyte  respectively.  The  values  of  R  were  estimated  from  the 
complex  impedance  plots  of  the  symmetric  cell-800.  The  Arrhenius 
plot  for  GDC  is  shown  in  Fig.  13.  The  activation  energy  Ea  =  0.88  eV 


T 
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Fig.  8.  Variation  of  log  or  at  600  °C  and  activation  enthalpy  (Ea)  with  x  for 
Nd2_xCexCu04±5. 
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Fig.  9.  Variation  of  d.c.  conductivity  at  650  °C  for  Ndi.8Ce0.2CuO4±5  with  Po2- 
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Fig.  10.  Complex  impedance  plots  for  cell-800  at  (a)  660  and  (b)  480  °C. 
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Fig.  11.  Complex  impedance  plots  at  various  temperatures  for  cell-1200. 
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Fig.  13.  Arrhenius  plots  for  bulk  conductivity  of  GDC  electrolyte  obtained  from  EIS 
on  cell-800. 


determined  in  the  present  work  is  in  good  agreement  with  the 
earlier  reported  value  of  Ea  =  0.9  eV  [28]. 

A  comparison  of  the  temperature-dependent  ASR  of  cell-800 
and  cell-1200  is  given  in  Fig.  14.  A  transition  at  approximately 
620  °C  is  clearly  seen  in  both  cells.  The  observed  transition  tem¬ 
perature  (Tc)  matches  closely  with  the  Tc&  625  °C,  corresponding 
to  NTC  to  PTC  transition  (Fig.  7).  The  ASR  below  and  above  Tc  obeys 
the  Arrhenius  law.  Furthermore,  the  ASR  and  activation  energy  in 
the  entire  temperature  range  of  the  measurement  are  higher  for 
cell-1200  than  for  cell-800.  The  sintering  of  the  cell  at  a  higher 
temperature  (1200  °C)  gives  rise  to  grain  growth,  thereby  reduc¬ 
ing  the  number  of  pores  in  the  cathode  (Fig.  4).  Both  these  factors 
reduce  the  catalytic  activity  and  result  in  the  high  ASR  and  activa¬ 
tion  energy  for  cell-1200.  The  activation  energies  in  the  low-and 
high-temperature  regions  are  compared  in  Table  2.  The  transition 
in  the  temperature-dependent  ASR  at  550  °C  has  been  observed 


Fig.  12.  Electrical  equivalent  circuit  for  symmetric  cells. 


Fig.  14.  Arrhenius  plots  of  ASR  of  cell-800  and  cell-1200. 
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Table  2 

Comparison  of  ASR  values  determined  in  this  study  with  values  reported  in  the 
literature. 
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Sr# 

Cathode 

ASR  (ohm cm2) 

T(°C) 

Ea  (eV) 

Reference 

1 

NduSro.3  CUO4 

0.24 

700 

[11] 

2 

La2Nio.6Cuo.404+($ 

1.0 

850 

[19] 

3 

Ndi.8Ceo.2Cu04 

- 

- 

1.06(H) 
1.44  (HT) 

[10] 

4 

Ndi.8Ceo.2Cu04 

0.13 

700 

0.95  (LT) 
1.73  (HT) 

a 

LT  -  low  temperature  region;  HT  -  high-temperature  region. 
a  Present  study 


by  Soorie  and  Skinner  [10].  The  disparity  between  the  values  of  Tc 
and  Ea  of  the  present  study  and  reported  values  [10]  (Table  2)  is 
attributed  to  the  variation  in  grain  size  and  grain  morphology  that 
result  from  differences  in  the  preparative  methods. 

The  variation  of  the  peak  frequency  (fp)  corresponding  to  the 
low-frequency  semicircle  in  the  complex  impedance  plane  with 
temperature  is  depicted  in  Fig.  15  and  suggests  Arrhenius-like 
behavior.  The  change  in  slope  at  620  °C  coincides  with  the  tran¬ 
sition  observed  in  the  temperature-dependent  ASR.  The  activation 
energy  (Ea)  is  determined  using  the  following  relationship: 

/p=/0exp(^),  (9) 

where  /p,  and  f0  are  the  jump  frequency  and  attempt  frequency 
respectively.  The  Ea  =  0.95  and  1 .73  eV  which  correspond  to  the  low- 
and  high-temperature  regions,  respectively,  and  obtained  from  the 
relaxation  frequency  plot  are  in  good  agreement  with  those  deter¬ 
mined  from  the  ASR  plot  (Table  2).  These  results  indicate  that  the 
charge  transfer  within  the  electrode  as  well  as  across  the  elec¬ 
trode/electrolyte  interface  took  place  via  the  hopping  mechanism. 

The  complex  impedance  plots  at  650  °C  for  cell-800  under  dif¬ 
ferent  d.c.  biasing  potentials  are  depicted  in  Fig.  16(a).  In  general, 
the  ASR  is  reduce  with  an  increased  biasing  potential  (Fig.  16(a) 
and  (b)).  On  the  other  hand,  there  is  a  negligible  change  in  the  elec¬ 
trolytic  bulk  resistance  with  the  variation  in  d.c.  biasing  potential, 
as  expected.  The  decrease  in  the  ASR  with  increased  d.c.  potential 
is  in  agreement  with  previous  studies  [29]. 

The  complex  impedance  plots  for  different  values  of  Pq2  for  cell- 
800  are  shown  in  Fig.  17(a).  The  x-axis  intercept  corresponding  to 
the  low-frequency  semicircular  arc  changes  with  Pq2.  On  the  other 
hand,  the  electrolytic  bulk  resistance/conductivity  is  least  affected 
by  changes  in  Pq2.  The  linear  dependence  of  log(ASR)  on  log(Po2) 


Z'  /  ohm  — ► 


DC  potential  /  V — ► 

Fig.  16.  (a)  Complex  impedance  plots  of  cell-800  at  650  °C  under  different  d.c.  polar¬ 
ization  potential  and  (b)  variation  of  ASR  with  d.c.  polarization  potential. 
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Fig.  15.  Arrhenius  plots  of  electrode  relaxation  frequency, /p,  obtained  from  cell- 
800. 


Fig.  17.  (a)  Complex  impedance  plots  at  650  °C  at  different  oxygen  partial  pressures 
and  (b)  variation  of  ASR  with  Pq2  for  cell-800. 
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shown  in  Fig.  17(b)  can  be  expressed  as: 

ASR  =  ASR0  (Pq2  )n  (10) 

The  value  of  n  gives  useful  information  about  the  type  of  species 
involved  in  the  reactions  at  the  electrode  [30]. 

n  =  1,  02(g)^02,abs  (11) 

n  —  2’  ^2,abs  ^  20abs  (12) 

n=l  Oabs  +  2e'  +  V0"^Oox  (13) 


The  value  of  n  =  0.23  estimated  in  the  present  study  indicates 
the  absorption  of  oxygen  in  the  oxide  lattice,  giving  rise  to  the  for¬ 
mation  of  electrons  and  oxygen  vacancies  according  to  Eq.  (13). 
The  increase  in  concentration  of  such  oxygen  vacancies  and  elec¬ 
trons  in  the  oxide  lattice  with  increased  Po2  reduces  the  ASR 
(Fig.  17(b)). 

The  ASRs  obtained  in  the  present  study  is  compared  with 
the  reported  values  in  Table  2.  A  close  scrutiny  of  Table  2 
reveals  that  the  Ndi.8Ce0.2CuO4±5  cathode  exhibits  the  low¬ 
est  ASR  compared  with  the  cathodes  reported  earlier.  The  low 
ASR  for  Ndi.8Ceo.2Cu04±5  in  the  present  study  is  due  to  its 
nanocrystalline  nature.  In  fact,  the  nanocrystalline  material  not 
only  enhances  the  electrical  conductivity,  but  also  improves  the 
electrode/electrolyte  interface  contact.  Additionally,  the  surface- 
to-volume  ratio  increases  and  leads  to  increased  surface  area  for 
catalytic  activity.  All  these  factors  are  responsible  for  the  reduction 
in  the  ASR. 

4.  Conclusions 

The  combustion  technique  considerably  reduces  the  sin¬ 
tering  temperature  and  time  and  also  yields  nanostructured 
Nd2_xCexCu04±($  solid  solutions.  Solid  solubility  up  to  x  =  0.2 
is  established  by  X-ray  powder  diffraction.  The  temperature- 
dependent  conductivity  exhibits  a  PTC  to  NTC  transition  at 
625  °C  for  Nd2_xCexCu04±(5  within  the  solid  solubility  limit.  The 
Ndli8Ce0.2QiO4±(5  cathode  is  thermochemically  stable  against  GDC, 
though  it  reacts  chemically  with  YSZ  at  high  temperature.  The 
catalytic  activity  of  Ndi.8Ceo.2Cu04±($  is  established  from  an  EIS 
study.  Sintering  of  the  symmetric  cell  at  high  temperature  increases 
the  ASR  due  to  grain  growth.  Ndi.8Ce0.2CuO4±($,  prepared  using 
the  combustion  technique,  promises  to  be  a  potential  cathode  for 
ITSOFC  applications. 


Acknowledgments 

The  authors  are  thankful  to  UGC,  New  Delhi  and  DST,  New 
Delhi  for  financial  support  through  the  SAP/DRS  and  FIST  programs, 
respectively,  to  carry  out  this  work.  Ms.  A.  P.  Khandale  is  grate¬ 
ful  to  UGC,  New  Delhi  for  the  award  of  a  Rajiv  Gandhi  fellowship 
(F14-2(SC)/2007(SA-III)). 

References 

[1  ]  M.  Dokiya,  Solid  State  Ionics  152/153  (2002)  383-392. 

[2]  C.  Chervin,  R.S.  Glass,  S.M.  Kauzlarich,  Solid  State  Ionics  176  (2005)  17-23. 

[3]  H.  Ullmann,  N.  Trofimenko,  F.  Tietz,  D.  Stover,  A.  Ahmad-Khanlou,  Solid  State 
Ionics  138  (2000)  79-90. 

[4]  X.Y.  Xu,  C.R.  Xia,  G.L.  Xiao,  D.K.  Peng,  Solid  State  Ionics  176  (2005)  1513-1520. 

[5]  V.V.  Kharton,  A.P.  Viskup,  A.V.  Kovalesvsky,  E.N.  Naumovich,  F.M.B.  Marques, 
Solid  State  Ionics  143  (2001 )  337-353. 

[6]  X.  Dong,  Z.  Wu,  X.  Chang,  W.  Jin,  N.  Xu,  Ind.  Eng.  Chem.  Res.  46  (2007) 
6910-6915. 

[7]  D.J.  Buttrey,  P.  Ganguly,  J.M.  Honig,  C.N.R.  Rao,  R.R.  Schartman,  G.N.  Subbanna, 
J.  Solid  State  Chem.  74  (1988)  233-238. 

[8]  J.D.  Jorgensen,  B.  Dabrowski,  S.  Pei,  D.R.  Richards,  D.G.  Hinks,  Phys.  Rev.  B  40 
(1988)2187-2199. 

[9]  A.  Mehta,  P.J.  Heaney,  Phys.  Rev.  B  49  (1994)  563-571. 

[10]  M.  Soorie,  S.J.  Skinner,  Solid  State  Ionics  177  (2006)  2081-2086. 

[  1 1  ]  X.  Ding,  X.  Kong,  J.  Jiang,  C.  Cui,  J.  Hydrogen  Energy  34  (2009)  6869-6875. 

[12]  A.  Aguadero.J.A.  Alonso,  M.T.  Fernandez,  M.J.  Escudero,  L.  Daza,J.  Power  Sources 
169(2007)  17-21. 

[13]  J.  Wan,  J.B.  Goodenough,  J.H.  Zhu,  Solid  State  Ionics  178  (2007)  281-286. 

[14]  M.J.  Escudero,  A.  Aguadero.J.  A.  Alonso,  L.  Daza,J.  Electroanal.  Chem.  611  (2007) 
107-116. 

[15]  M.A.  Daraukh,  V.V.  Vashook,  H.  Ullmann,  F.  Tietz,  I.  Arual  Raj,  Solid  State  Ionics 
158(2003)  141-150. 

[16]  Q.  Li,  H.  Zhao,  L.H.  Huo,  L.  Sun,  X.L.  Cheng,  J.C.  Grenier,  Electrochem.  Commun. 
9(2007)  1508-1512. 

[17]  V.V.  Kharton,  A. A.  Yaremchenko,  E.N.  Naumovich,  J.  Solid  State  Electrochem.  3 
(1999)303-326. 

[18]  S.F.  Palguev,  V.K.  Gilderman,  V.I.  Zemstov,  High  Temperature  Oxide  Electronic 
Conductors  For  Electrochemical  Devices,  Nauka,  Mascow,  1990. 

[19]  A.  Aguadero,  J.A.  Alonso,  M.J.  Escudero,  L.  Daza,  Solid  State  Ionics  179  (11-12) 
(2008)393-400. 

[20]  V.I.  Voronin,  A.E.  Kar’kin,  B.N.  Goshchitski,  Phys.  Solid  State  40  (2)  (1998) 
157-162. 

[21]  J.B.  Goodenough,  A.  Manthiram.J.  Solid  State  Chem.  88  (1990)  115-139. 

[22]  P.  Ganguly,  C.N.R.  Rao,  J.  Solid  State  Chem.  53  (1984)  193-216. 

[23]  G.  Spinolo,  M.  Scavini,  P.  Ghigna,  G.  Chiodelli,  G.  Flor,  Physica  C  254  (1995) 
359-369. 

[24]  T.J.B.  Holland,  S.A.T.  Redfern,  Miner.  Mag.  61  (1997)  65-77. 

[25]  J.S.  Kim,  E.P.  Kvam,  Physica  C  292  (1997)  203-210. 

[26]  H.  Takagi,  S.  Uchida,  Y.  Tokura,  Phys.  Rev.  Lett.  62  (1989)  1197-1199. 

[27]  O.G.  Singh,  B.D.  Padalia,  O.  Prakash,  K.  Suba,  A.V.  Narlikar,  L.C.  Gupta,  Physica  C 
219(1994)156-161. 

[28]  K.  Singh,  S.A.  Acharya,  S.S.  Bhoga,  Ionics  12  (2006)  295-301. 

[29]  F.  Mauvy,  J.M.  Bassat,  E.  Boehm,  J.P.  Manaud,  P.  Dordor,  J.C.  Grenier,  Solid  State 
Ionics  158(2003)  17-28. 

[30]  E.  Siebert,  A.  Hammouche,  M.  Kleitz,  Electrochim.  Acta  40  (1995)  1741-1753. 


